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We use the moment approach of Nolting (exact sum rules) (Z. Physik 255, 25 (1972)) for the 
attractive Hubbard model in the superconducting phase. Our diagonal and off - diagonal spectral 
' functions are constructed and evaluated with the sum rules. They reduce to the BCS limit for 

Q\ , weak interaction. However, the presence of correlations modify the BCS picture dramatically. For 

example, due to the presence of correlations we have postulated a three - pole ansatz for the diagonal 
Green function, G(k,ui), while the off - diagonal one, B(k,cu), is supposed to have two poles. In 
the paper we present results for the three spectral weights of the diagonal Green function, ay(fc), j 
= 1,2,3. Our results compare reasonably well with more elaborated auto - consistent highly non - 
linear equations (double fluctuation calculations in the T— Matrix approach of one of the authors). 
■ Then, the physical picture which emerges is that the lower Hubbard band is split due to the super- 

conducting gap and the upper Hubbard band remains mostly unmodified. 
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I. INTRODUCTION 



The study of correlations has been renewed again by the the discovery of high T c superconducting oxide materials 
(HTSC), jl| since that these materials exhibit a short coherence lenght, £, and a very large penetration depth, A, 
which become them almost extreme type-II superconductors, i.e. ke|>1. 

In order to study the effect of correlations we will use the on-site attractive Hubbard Hamiltonian put forward by 



Micnas et al || as a phenomenological model for describingthe HTSC. Previous authors have used this model to 
study the bismuthate superconductors. H Denteneer et al [Q recognize that if the phase diagram of the Hubbard 
model were fully understood, it might form the basis of understanding correlated electrons as much as the Ising 
model did for understanding critical phenomena. More recently, Schneider et al || have applied this model to explain 
universal properties of several families of these compounds, like the relation between transition temperature, magnetic 
' penetration depth and gap, at zero temperature. 

In this paper, we will use the exact relations of Nolting H for the one-particle diagonal spectral function together 
with the exact relations for the one-particle off-diagonal spectral function, i.e., the anomalous Green's function, to 
study systems with broken symmetry. ffo 
The model we study is the following 

H = -t 2J Aa c i" + u ^2 n it n H ~ A*zJ nta ' W 

<i,j>(7 i ia 

where cl a (ci a ) are creation (annihilation) electron operators with spin a. rii a = cl a Ci a , t is a hopping matrix element 
between the nearest sites (n.n.) i and j, and U is the onsite interaction. \i is the chemical potential and we work in 
the grand canonical ensemble. In the present study we consider an attractive interaction, U = — | U | < 0. The 
Hamiltonian of Eq. ([j]) has been studied in detail by Micnas et al ||. Nolting || and Kalashnikov and Fradkin j?]] 
have derived exact relations for the one-particle spectral functions which are given in Ref. |j. 
We postulate for the diagonal spectral function, A(k, u>), the following ansatz 

A{k,tu) = ai (k)S{uj - E % ) + a 2 {k)5{io + E % ) + a 3 {k)S(uj - %) (2) 

and the off-diagonal spectral function, B(k,u>), is supposed to have two peaks (BCS-like), as follows 
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B(k, u;) = y/ ai (k)a 2 (k) [5{lo - E % ) - 5(lu + E % )] , Q^k) = E % , 2 (fc) = -E % 



(3) 



Our ansatz is based on the assumption that the diagonal spectral function is modified by the presence of correlations 
(we define correlations when we have the regime U/W ~ 1, where W — 8t is the bandwidth in two dimensions). This 
assumption has been numerically checked to be correct in Ref. |9| where the authors have obtained three bands for 
the diagonal spectral function and two bands for the off-diagonal spectral function after performing the analytical 
continuation. In Eq. (Q), the first two poles represent the behavior around the chemical potential and the third pole 
is due to the influence of the upper Hubbard band, which in the case of the atomic limit is the upper Hubbard band 
or the band of single occupied states. The lower band, now split in two, is the band of double occupied states for the 
case of allmost atomic limit. Please, see reference M], where the attractive Hubbard model has been studied in the 
normal phase with two different methods: T-Matrix and moments. 

The order parameter, A (T), and the chemical potential (fixed electron density), p, are calculated from the follow- 
ing equations: 



fc fc 

For example, the quasi-particle energy, E^, has a gap. As it was discussed in Ref. the gap is expressed in terms 
of the order parameter as 

A(r,fc) = A (T) 2E ' k , (5) 
Ej: + ai(k) 

concluding from Eq. (^J) that our gap equation, which is a manifestation of singularities in the density of states, is 
fc-dependent. Going back to our local Hamiltonian (see Eq. (|l|)), which in reciprocal space is a constant, then at the 
mean field level we should obtain a pure s-wawe, while going beyond mean field approximation we have been able to 
derive a fc-dependent gap. In our case, by including correlations, i.e., the selection rules for a 2 and 03 (See Ref. || for 
more details) we have been able to modify the BCS results, from pure s-wave to a wave vector dependent energy gap. 
Let us point out that in the attractive Hubbard model we always find s— type wave symmetry as it has conclusively 
shown in Ref. in the Hubbard-I approximation for both the diagonal and the off-diagonal one-particle Green 
functions. 

We would like to call the reader's attention that our Ansatz is based on the assumption that the role of correlations 
is mainly taken into account in the diagonal one-particle spectral function. This implies that the off-diagonal order 
parameter has been taken at the BCS level, i.e., to M^i/k). In a previous work, |l2| it is found that both the 
diagonal and off-diagonal spectral functions have four peaks, symmetric in pairs, for U/t = —4. However, when we 
have included double fluctuations, i.e., fluctuations in both the diagonal and off-diagonal self-energies as in Ref. ||, 
we have three resolved peaks for A(k, ui) and two well resolved peaks for B(k, lu). So, Ref. Q has motived us to adopt 
three and two peaks for A(k, uj) and B(k, oj), respectively, for U/t = —4.0. We have restricted ourselves to wavevectors 
close to the chemical potential. If we want to consider two other peaks for the off-diagonal spectral function we have 
to include more moments, which is beyond the present work. Also, we have neglected life-time effects in Eqs. (|2 -p[) , 
which would require to approach the phase transition with damping in the self-energy as it has been done in Ref. [ 13) . 

The equations we have to solve are the following 

1 1 tanh(^l) 

U N ^ 2E k 1 ' 

k 

p = 1 J2 U(k)/(E k ) + a 2 (k)/(-£ k ) + ai f< \~ ^ W /(flh) l (7) 



In Fig. 1 we show the dependence of the diagonal spectral weigths, ctj(k), j=l,2,3 along the diagonal of the Brillouin 
zone. We can appreciate that the area of the spectral function is equal to one (ai + a 2 + a% = 1) by construction. We 
have taken U/t = —4.0, p = 0.1 and T/t = 0.01. After self-consistency of our highly non-linear equations we obtain 
p/t = —4.19, B/t = 1.664 and Ao(T) = 0.52. It is worth mention that p < —At which indicates that we are in the low 
density regime (Bose - Einstein limit). In Fig. 2 we have the diagonal spectral weights along the fc^-axis (k = (fc,0)). 
What we observe is a strong variation of the ctj(k) with direction. In Fig. 3 we have the fc-dependence of 7j(fc), j = 
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1,2 along the two directions as previously discussed in Figures 1 and 2. These -fj(k)'s are defined by Eq. (19) in Ref. 
||. They are highly ^-dependent. The meaning of these variables is that they represent corrections to the diagonal 
sum rules (ei2(fe) and 0,3 (fc)). As we have said in Ref. S these corrections are important, i.e., of the order of 20 - 25 
higher corrections on the order parameter (or equivalently, 7j(fc)'s). For T T c these corrections are not important, 
though. 

By using the moment approach in the presence of off-diagonal long range order, i.e.,A Q (T) ^ 0, to the negative 
Hubbard model, we have worked out the superconducting phase calculated with the use of moments (exact sum 
rules for both the diagonal and off-diagonal spectral densities) with three and two and two peaks, respectively. The 
physical meaning of each of these peaks has been discussed. This assumption is equivalent to solve the diagonal 
spectral weights in powers of the order parameter square. At the same time, we have seen that the effect of the third 
band is to renormalize the order parameter producing an energy gap which is fc-dependent. In consequence, in our 
approach the order parameter, A D (T), and the energy gap, A(T, k) are not the same. We would like to end by saying 
that there are strictly mean - field calculations (even if they go beyond Hartree - Fock) as the ones done by Martin 
- Rodero and Flores JTJ] and Kuchiev and Sushkov [jlB). Our approach presented here is able to capture the double 
fluctuaction calculation of Ref. This clearly indicates that we are beyond a simple BCS calculations. 

One of the authors ( J JRN) would like to acknowledge partial support from CONDES-LUZ and also from CONICIT 
(project F-139). We thank Maria Dolores Garcia for a reading of the manuscript. 
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Figures 



Figure 1. The spectral weights, ctj(k), j = 1, 2, 3 along the diagonal of the Brillouin zone. The first diagonal sum rule 
is satisfied. U/t = -4.0, p = 0.1 and T/t = 0.01. 

Figure 2. The spectral weights, Uj(k), j = 1, 2, 3 along the k x direction of the Brillouin zone. Same parameters as in 
Figure 1. 

Figure 3. 7j(fc), j = 1, 2, 3 along the two directions discussed in Figures 1 and 2. Same parameters as before. 
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°— ctjOc) |x/t = -4.19 V 1 = °- 52 
A _ a (k) B/t = 1 .664 u/t = -4.0 
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